Nonequilibrium transition induced by mass media in a model for social influence 
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We study the effect of mass media, modeled as an applied external field, on a social system based 
on Axelrod's model for the dissemination of culture. The numerical simulations show that the sys- 
tem undergoes a nonequilibrium phase transition between an ordered phase (homogeneous culture) 
specified by the mass media and a disordered (culturally fragmented) one. The critical boundary 
separating these phases is calculated on the parameter space of the system, given by the intensity of 
the mass media influence and the number of options per cultural attribute. Counterintuitively, mass 
media can induce cultural diversity when its intensity is above some threshold value. The nature of 
the phase transition changes from continuous to discontinuous at some critical value of the number 
of options. 

PACS numbers: 89.75.Fb, 87.23. Ge, 05.50.+q 



In recent years, the research on complex systems has 
extended to social science in order to understand how 
collective behaviors arise in social systems. Several 
mathematical models, inspired by analogies with phys- 
ical systems, have been proposed to describe a variety 
of phenomena occurring in social dynamics 0, 0, H, H|- 
Processes such as self-organization, cooperation, epi- 
demic spreading, opinion formation, propagation of in- 
formation, economic exchanges and evolution of social 
structures have been studied by means of discrete-time, 
discrete-space dynamical systems. In this context, there 
has been interest in the model introduced by Axelrod |f| 
to investigate the dissemination of culture among inter- 
acting agents in a social system 0, S S Eli El ■ From 
the point of view of statistical physics, this model is ap- 
pealing because it exhibits nontrivial out of equilibrium 
dynamics, as in other well studied systems with phase or- 
dering properties ^3 ■ Studies on this model have mainly 
focused on the collective properties that result from the 
interactions between the elements representing endoge- 
nous social influences. 

In this paper we investigate the effect of external cul- 
tural influences such as mass media on a social system. 
Our approach is based on the adaptive interaction dy- 
namics of Axelrod's model for the dissemination of cul- 
ture. Agents can interact with their neighbors in the sys- 
tem and with the mass media according to the cultural 
similarities that they share, in each case. The concept of 
culture is intended here as a set of individual features or 
attributes that are subject to social or external influence. 
The numerical simulations show that, depending on the 
value of a parameter that represents the intensity of the 
mass media influence and on the number of options avail- 
able per cultural attribute, the system displays a phase 
transition between an specific ordered phase (a homoge- 
neous culture) imposed by the mass media and a disor- 
dered (culturally fragmented) phase. Surprisingly, mass 
media can induce cultural diversity when its intensity is 
above some threshold value. The nature of this transi- 
tion changes from continuous to discontinuous when the 



number of options per cultural feature is increased. 

The model consists of N agents as the sites of a 
square lattice. The cultural state Cj of agent i is de- 
fined as a vector of F components (cultural features) 
Ci = (cjij Cj2, ■ ■ ■ , oif)- Each a*/ can take any of the q 
values in the set {0, 1, . . . , q— 1} (cultural traits), initially 
assigned randomly with a uniform distribution. There 
are q F possible cultural vectors. We define a mass me- 
dia cultural message as a vector M = (/ii,/i2, ■ ■ ■ iH-f), 
where \if G {0, 1, . . . ,q — 1}, that can interact with any 
of the agents in the system. We also define a parame- 
ter B G [0, 1] that measures the relative intensity of the 
mass media message with respect to the local interac- 
tions, or the probability that this message has to attract 
the attention of the agents in the system. The param- 
eter B represents enhancing factors of the transmitted 
message that can be varied externally, such as its ampli- 
tude, frequency, attractiveness, etc. It is assumed that 
B is uniform, i.e., the mass media message reaches all 
the agents with the same intensity, as a uniform field. At 
any given time, we assume that any agent can either in- 
teract with the mass media message or with other agents 
in the system. Thus each agent in the network possesses 
a probability B of interacting with the message and a 
probability (1 — B) of interacting with its neighbors. 

The discrete-time dynamics of the system subject to 
the mass media influence is defined by iterating the fol- 
lowing steps: 

(1) Select at random an element i in the lattice (active 
element) having a cultural state c;. Its attention is drawn 
to interact with the message M with probability B. 

If the attention of element i is drawn to the message, 
then 

(2) Calculate the cultural overlap (number of shared 
features) between the active element and the message 

l(i,M) = J2f=i^if,f. r 

(3) If < l(i,M) < F, the element i and the mes- 
sage interact with probability l(i,M)/F. In case of in- 
teraction, choose h randomly such that <Jih ^ fJ-h and set 
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If the attention of element i is not caught by the mes- 
sage, then 

(4) Select at random a site j in the nearest neighbor- 
hood of site i. 

(5) Calculate the cultural overlap = 

(6) If < < F, sites i and j interact with proba- 
bility l(i,j)/F. In case of interaction, choose h randomly 
such that (Tih ^ (Tjh and set Oih = <Jjh- 

In this model, the fixed mass media cultural message 
M can interact with any element in the system with the 
same intensity or probability B, but the effect of that 
interaction may be different on each element, depend- 
ing on the specific cultural overlap between the element 
and the message. For simplicity, we are assuming a fixed 
mass media cultural message acting uniformly over the 
system, as for example in a global broadcasting; however 
several variations of this condition can be implemented 
by extending this basic algorithm. The case B = cor- 
responds to the original Axelrod's model. 

In any finite network the dynamics settles into an ab- 
sorbing or frozen state, characterized by either = 
or = F, Mi, j. Homogeneous or monocultural states 
correspond to = F, Vi, j, and obviously there are 
q F possible configurations of this state. Inhomogeneous 
or multicultural states consist of two or more homo- 
geneous domains interconnected by elements with zero 
overlap. A domain is a set of contiguous sites with iden- 
tical cultural traits. In the absence of external influences 
it has been shown that the system reaches ordered, mono- 
cultural states for q < q a and disordered, multicultural 
states for q > cu, where q Q is a critical value that de- 
pends on F 0, 0, lJ ■ This order-disorder phase tran- 
sition becomes better defined as N increases and it has 
been argued that it is of first order in two-dimensional 
systems 0, Q , while this transition is of second order in 
one-dimensional systems ^lj. The critical value of the 
number of traits for F — 10 has been numerically esti- 
mated at q a « 55 in two dimensions ||. 

When the mass media cultural influence is applied to 
the system, the order-disorder phase transition persists, 
but the critical value q for which the transition takes 
place decreases as the intensity of the message B is in- 
creased. Figure 1 shows the spatial configurations of the 
final states of the system subject to mass media, with 
F = 10 and q = 35 < q . In the absence of external 
message, i.e. B = 0, the system settles into any of the 
possible q F monocultural states. When the intensity of 
the message is increased, the system is driven towards a 
monocultural state equal to the state of the mass media, 
i.e., Ci = M, Vi. However, there is a critical value of the 
intensity B c ks 0.05 above which the system no longer 
converges to the state of the message M but reaches a 
multicultural state consisting of an increasing number of 
domains as B is increased. Domains having a cultural 
state equal to the mass media survive above the thresh- 
old intensity B c , but they become smaller in size as B is 
increased. Thus, we find the counterintuitive result that, 



above some threshold value of intensity, mass media ac- 
tually promotes cultural diversity in the system. 




FIG. 1: Spatial patterns for different values of the intensity 
B, for F = 10, q = 35, and N = 50 x 50. The color code of 
the mass media state M is white. Top left: B — 0; top right: 
B = 0.005; bottom left: B = 0.1; bottom right: B = 0.9. 

In order to characterize the transition from the mono- 
cultural state imposed by the message M to a multicul- 
tural state when the intensity B is varied, we consider 
as an order parameter the average fraction of cultural 
domains g = (N g )/N, where N g is the number of do- 
mains formed in the final state for a given realization of 
initial conditions [l^. Figure 2 shows g as a function of 
B for different values of the number of traits q. For val- 
ues of q > q , the system always reaches a multicultural 
state (with (N g ) 3> 1), independently of the intensity 
B of the message. On the other hand, for each value of 
q < q , we observe a transition at a critical value B c , from 
the homogeneous cultural state imposed by the message, 
characterized by values g <tC 1, to a multicultural state 
for which g increases with B. The critical value B c de- 
creases with increasing q for q < q Q . The variation of the 
order parameter g near the value B c can be character- 
ized by a critical exponent (3(q) as the scaling relation 
g ~ (B — Bc) 13 ^' . Figure 2 also shows log-log plots of g 
vs. (B — B c ) for different values of q < q . The exponent 
P as a function of q can be calculated from the slope of 
each curve. 

In Fig. 3 we show the resulting graph of /3 vs. q, for 
q < q corresponding to F = 10. The dependence of the 
exponent f3 with q is well accounted by the linear relation 
P{q) (x (q — q). As q increases towards the value q 0l 
the exponent /3 becomes smaller and the corresponding 
phase transition from the monocultural state induced by 
the message M to multicultural state gets more abrupt. 
The character of the phase transition changes from higher 
order (/? > 1 ) to second order (0 < j3 < 1) at q m 30. 
At the value q = q for which the exponent (3 vanishes, 
the character of the transition should change from second 
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FIG. 2: Order parameter j as a function of the intensity B 
for different values of q, with F = 10. Error bars indicate 
typical standard deviations about mean values obtained over 
10 realizations for each value of B. The inset shows a log- 
log plot of g vs. (B — B c ) for different values of q. Size 
TV = 50 x 50. Values of q are 20 (squares), 30 (diamonds), 45 
(circles), and 60 (triangles). 



order to first order. Figure 3 shows the extrapolation of 
the straight line (3{q) until its intersection with the q axis, 
predicting a critical value q = 54 for the occurrence of 
a first order transition. 
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FIG. 3: Critical exponent f3 for the phase transition induced 
by mass media, as a function of q. 

The critical values of the intensity of the message B c 
as a function of q, calculated from Fig. 2, are shown in 
Fig. 4. The points in Fig. 4 arc well fitted by the re- 



lation B c oc (e 
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), which yields B c = for 



q = q a = 54 as indicated in that figure. Thus, in ab- 
sence of mass media, the critical value of the number of 
traits for the first order phase transition from cultural 
homogeneity to a multicultural state predicted by our 
model is q Q = 54. This agrees with the approximated 
critical value obtained numerically in Ref. Q for the 



order-disorder transition in the original Axclrod model 
with F = 10. The critical curve B c vs. q in Fig. 4 sepa- 
rates the multicultural region from the region where the 
monocultural state induced by the mass media occurs on 
the space of parameters (B,q). This critical curve re- 
flects the competition between the tendency towards the 
specific ordering imposed by the mass media and the dis- 
order present in the initial configuration of the system, 
represented by the number of options per feature q. Fig- 
ure 4 reveals that the ability of the mass media message 
to transmit its state on the entire system decreases as the 
number of options per cultural feature increases. Above 
this threshold curve, the mass media is not longer able 
to drive the system towards homogeneity, and a multi- 
cultural state sets in. The above results suggest that the 
order parameter g behaves, near the critical boundary, 



as g(B, q) ~ [B - k 2 (e 
ki and k% are constants. 



where k\ , 
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FIG. 4: Critical boundary B c vs. q separating the region 
of the monocultural phase induced by mass media from that 
of multicultural states, as indicated. The dashed line corre- 
sponds to the fitting of the points. Fixed F = 10. 

In the presence of the mass media message, the system 
reaches frozen configurations of each phase (homogeneous 
or multicultural) much faster than in the original Axclrod 
model having only local interactions. In Fig. 5 we show 
a log-log plot of the average time r to reach a frozen 
configuration as a function of the intensity B, for several 
values of q < q a . For B < B c , the time r decays accord- 
ing o r ~ B~ 95 , independently of q for q < q a (where 
the exponent is determined from the slopes of the curves 
in Fig. 5). The times for convergence to the monocul- 
tural state induced by the message are much longer than 
the corresponding times to reach a multicultural state 
beyond the critical boundary in Fig. 4. Therefore, for 
values of intensity B < B c the mass media message is 
able to transmit its cultural state to all the elements in 
the system, but the process takes longer times than the 
partial convergence to the message state achieved when 
B > B c . 
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FIG. 5: Log-log plot of the average transient time to freezing 
r as a function of B, for different values of q < q . Here 
q — 10 (pentagons), q = 20 (squares), q — 30 (diamonds). 

The interaction with the mass media message M accel- 
erates the convergence of those vectors Cj in the system 
that possess an overlap < l(i,M) < F towards the 
cultural state prescribed by M. The dynamics of the 
interaction with the mass media is essentially a coarsen- 
ing process of such domains. For values of B and q be- 
low the critical boundary, this process takes longer times 
and allows enough local interactions to take place and 
spread the cultural state of the message, leading to the 
growth and coalescence of domains with the state M into 
a single domain of size comparable to that of the whole 
system. For values (B, q) above the critical boundary, 
the process of convergence to the message and the for- 
mation of domains with this cultural vector occurs faster; 
this contributes to the early distinction of the domains 
with state M from the other domains and to limit their 
growth. Thus, the rapid convergence towards the mass 
media message when its intensity is above a threshold 
value ends up inducing cultural diversity in the system. 
The transition between the two phases in Fig. 4 can also 



be seen as separating two different dynamical regimes; 
one having a fast decay time for values (B,q) above the 
critical boundary and another characterized by a slow 
decay for values (B, q) below that boundary. 

In summary, we have presented a model for the influ- 
ence of mass media on a social system as an extension 
of Axelrod's model for the dissemination of culture. The 
mass media cultural message has been assumed as a fixed 
vector acting uniformly over the system. We have found 
the nontrivial result that mass media can actually induce 
cultural diversity in conditions for which it was absent, 
as for values q < q Q in the original Axelrod's model. We 
have verified that this effect persists in lattices with pe- 
riodic boundary conditions or with different local geome- 
tries. 

We have calculated the critical boundary in the space 
of parameters (B, q) that separates the monocultural 
phase induced by the mass media message from the multi- 
cultural phase. These two phases correspond to different 
dynamical regimes of the system as it evolves towards its 
frozen final state; the first having a slow decay and the 
second characterized by a fast decay. The character of 
the phase transition changes progressively from continu- 
ous to discontinuous as the number of traits per feature 
is increased. The critical value q for the onset of the 
first order phase transition was predicted from the linear 
dependence observed in the critical exponent /3 vs. q, as 
well as from the curve fitting the critical boundary. The 
predicted value q agrees with the value estimated previ- 
ously by numerical simulations in the original Axelrod's 
model. The scaling behavior of the order parameter g 
characterizing the transition near the critical boundary 
was also found numerically. Future extensions of this ba- 
sic model should include the consideration of mass media 
messages varying in time and/or space, competing mes- 
sages, noise, and more complex networks of interactions. 
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